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Abstract: The PSI low-energy /iSR spectrometer is an instrument dedicated to muon spin ro¬ 
tation and relaxation measurements. Knowledge of the muon beam parameters such as spatial, 
kinetic energy and arrival-time distributions at the sample position are important ingredients to an¬ 
alyze the /rSR spectra. We present here the measured energy losses in the thin carbon foil of the 
muon start detector deduced from time-of-flight measurements. Muonium formation in the thin car¬ 
bon foil (10 nm thickness) of the muon start detector also affect the measurable decay asymmetry 
and therefore need to be accounted for. Muonium formation and energy losses in the start detec¬ 
tor, whose relevance increase with decreasing muon implantation energy (< 10 keV), have been 
implemented in Geant4 Monte Carlo simulation to reproduce the measured time-of-flight spectra. 
Simulated and measured time-of-flight and beam spot agrees only if a small fraction of so called 
“unmoderated” muons which contaminate the mono-energetic muon beam of the /rSR spectrometer 
is introduced. Moreover the sensitivity of the beam size and related upstream-downstream asym¬ 
metry for a specially shaped “nose” sample plate has been studied for various beam line settings, 
which is of relevance for the study of thermal muonium emission into vacuum from mesoporous 
silica at cryogenic temperatures. 
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1. Introduction 

Polarized positive muons ean be used to investigate struetural properties and dynamieal pro- 
eesses of solid states via so ealled /rSR teehnique [1] which stands for a collection of methods as 
Muon Spin Rotation, Relaxation and Resonance. A low-energy beam with tunable energy in 
the keV regime [2] plays an important role in the /rSR field because these muons can be used to in¬ 
vestigate thin films. The Low Energy Muon (LEM) beam line at PSI shown in figure 1 is delivering 
/r+ with tunable kinetic energy between 0.5 and 30 keV allowing the study of thin films and multi¬ 
layers as a function of the implantation depth ranging from 0.5 nm to a few hundred nm [3, 4]. 
The low-energy are obtained by moderating a surface beam (4 MeV energy) from the /xE4 
beam line [5] with a 125 /xm thick Ag foil coated with a 200-300 nm thick layer of solid Ar-N 2 [6]. 
The moderated leaving the solid Ar-N 2 surface have mean kinetic energy of 15 eV. As the mod¬ 
erator is placed at high voltage (typically Vmod = 15 kV), after leaving the moderator the jj.^ are 
accelerated to about 15 keV kinetic energy. Using various electrostatic elements, the /r+ are trans¬ 
ported from the moderator to the sample region. Neglecting the various focusing elements (einzel 
lenses and conical lenses), the are first bent by a 45° electrostatic mirror and then transported 
through a spin rotator. Subsequently, they cross a thin carbon foil (C-foil) which acts as a start 
detector before being implanted into the sample. The nominal density and thickness of the C-foil 
are ~2 /rg/cm^ and 10 nm, respectively. 

While crossing the C-foil the is ejecting several electrons (on average 3) of few eV energy 
from the foils surface whose detection define the implantation time in the sample and the 
start of the event in the data acquisition. The initially mono-energetic arriving at the start 
detector undergoes energy and angular straggling which degrade the beam quality and affect the 
measured /rSR time spectra. Moreover, by traversing the thin C-foil, a fraction of the /d+ can 
undergo charge-exchange and leave the foil as muonium (Mu) or negatively charged muonium 
(Mu^) which decrease the measurable total muon decay asymmetry. 

Another complication is represented by the low-energy tails of the leaving the Ar-N 2 mod¬ 
erator not as epithermal at the eV energy but as only partially moderated muons whose energy 
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about 0.01% to an energy of about 15 eV before being re-accelerated again to energies up to 20 keV. The 
beam is then bent by a 45° electrostatic mirror before going through the spin rotator and the start detector 
and arriving at the sample plate mounted on the cold finger of the cryostat. 


is sufficiently low to be deflected by the 45° electrostatic mirror and transported to the sample 
region. We term these muons as “unmoderated” muons. 

All these processes, at the moderator and at the C-foil, affect the kinetic energy distribution of 
the leaving the C-foil and consequently the arrival time distribution at the sample position 
(relative to the signal in the start detector). This arrival time distribution needs to be known to 
understand the detailed shape of the decay positron time spectra at early times. 

In this paper, measurements of the energy loss in the thin C-foil at various energies which 
have been done using a time-of-flight (TOP) technique are presented. Energy losses and Mu for¬ 
mation in the 10 nm thick C-foil of the start detector as well as “unmoderated” fraction of /!+ 
have been implemented in the musrSim [8] simulation package which is based on Geant4 [9] to 
match the measured TOE spectra. The Geant4 simulation has been performed starting from the 
downstream of the moderator till the sample region. 

The EEM beam line was upgraded in 2012 to allow for longitudinal /iSR measurements, to 
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Figure 2. TOF spectrum from the start detector to a micro-channel plate (MCP2) detector placed at the 
sample plate position for 12 keV beam transport energy. The peaks are corresponding to the prompt 
photons, the C-foil Mu atoms and the respectively. The time axis is set with an arbitrary zero. The HV 
at the RA and L3 lenses are switched off. 

have a better suppression of the proton/ion background from the moderator and to have a better 
time resolution. The first two items have been achieved by the installation of a spin rotator [10] 
as shown in figure 1 after the electrostatic mirror, whereas the improved time resolution has been 
achieved by reducing the distance between the start detector and the sample region (from 1164 mm 
to 563 mm distance, c.f. figure 1 of [3]). 

The beam sizes at the sample position for variations of the beam line settings have been studied 
and compared with determinations obtained from the upstream-downstream asymmetry measure¬ 
ments. 

2. Energy loss in the thin C-foil via TOF measurements 

To determine the energy loss in the C-foil of the start detector a TOF technique was applied. A 
muon beam with well defined kinetic energy given by the moderator high voltage (HV) is focused 
into the C-foil of the start detector which is set to a negative HV (Vc = —3.38 kV). The secondary 
electrons knocked out in the downstream direction by the muon crossing the foil are first accelerated 
and then deflected by a system of grids towards a micro-channel-plate (MCP3) as shown in the 
inset of figure 1. The MCP3 signal provides the “start” signal of the TOF measurement. Another 
micro-channel-plate (MCP2) is placed at the position usually taken by the sample holder while 
performing the /fSR measurements. A signal from the MCP2 delivers the stop time of the TOF 
measurement. In addition, the detection of an e+ from muon decay in the positron counters of 
the /fSR spectrometer is required. The measured time spectra obtained in this way are shown in 
figure 2. Understanding these time spectra requires some knowledge of the processes occurring in 
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Figure 3. (a) Formation of Mu and Mu^ in the carbon foil due to the charge exchange process, (b) The 
charge state yields of the exiting the carbon foil as a function of incoming muon energy at the foil surface, 
according to a velocity scaling of the proton data parameterizations [11]. 


the C-foil. The muons can leave the foil at various charge states: Mu and Mu^, as depicted 

in figure 3(a). The equilibrium yield for these various charge states scaled from proton data [11] 
are shown in figure 3(b). As fhe model of H formafion (via overlap of fhe afomic slates and solid 
elecfron sfafes, and elecfron lunneling given in [11]) when a low-energy prolon beam crosses a 
Ihin C-foil depends only on fhe velocify, velocity scaling of fhe profon dala has been assumed fo 
calculale fhe Mu charge slate yield. 

The firsl peak al early limes (/peak) is ensued in Ihe following way: Ihe eleclrons knocked oul 
by a muon crossing Ihe C-foil are transported and detected in the MCP3 delivering the start time of 
the event. The avalanche process occurring in the channel walls of the MCP3 generates about 10^ 
electrons impinging with energies of a few hundred eV on the anode of MCP3, where they generate 
UV photons by ionization/recombination processes or Bremsstrahlung. Some of these photons 
could be detected (non-zero solid-angle acceptance and non-zero detection efficiency) in fhe MCP2 
providing Ihe slop lime of Ihe event The widlh of Ihe / peak is aboul 0.7 ns and il corresponds 
lo Ihe liming resolulion of Ihe TOP system (MCPs and eleclronics). ll was demonslraled lhal Ihe 
position of Ihis peak is nol affected by HV variations al Ihe RA conical lens or al Ihe L3 einzel lens 
or in Ihe slarl detector grids confirming lhal il is related wilh photon emission in Ihe MCP3. 

The second peak is caused by Mu atoms Iraveling from Ihe C-foil to Ihe MCP2 while Ihe Ihird 
(largesl) peak is caused by /i+. Also for Ihese peaks Ihe slarl signal of Ihe evenls is given by Ihe 
eleclrons emilled from Ihe C-foil. Since Ihe Mu motion is nol affected by Ihe eleclric fields in Ihe 
slarl detector, Ihe Mu peak position (relative to Ihe / peak) can be related in a simple way to Ihe 
kinetic energy of Ihe Mu after Ihe C-foil. The delayed timing of Ihe /i+ peak compared to Ihe Mu 
peak is mainly due to Ihe facl lhal /r+ has to overcome Ihe negative eleclroslalic potential in Ihe 
slarl detector. Thus, /r+ leave Ihe slarl detector al a smaller kinetic energy compared to Mu atoms 
producing Ihe observed delay. 

The time to al which a /i+ is crossing Ihe C-foil is given by: 

to = tY — AIfe — Ac , (2.1) 
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where ty is the position of the 7 peak, 

AtpE is the TOF of the knoeked out eleetrons from the C-foil to MCP3 and Ate = 1.67 ns is 
the TOF of a partiele with the speed of light from the MCP3 to the MCP2. A AtpE = 13.5(5) ns 
has been determined from the time differenee of two prompt peaks in the deeay time speetra in 
regular /rSR measurements. The earlier peak is eaused by positrons hitting the C-foil and produeing 
foil eleetrons sueeessively deteeted in the eounters surrounding the sample region and the later 
peak is eaused by positrons hitting the MCP3 direetly and then deteeted in the eounters. The 
TOF of foil eleetrons to the MCP3 is independent of the emission position at the C-foil but 
has a small dependenee on the emission angle. However, this effeet and the possible variation of 
the emission energy are ineluded in the uneertainty of AtpE. The TOF of through TD was 
simulated for various initial position on the C-foil. The resulting TOF spread is about 0.1 ns. Sinee 
the motion of the Mu atom, being a neutral system, is not affeeted by the eleetrie field in the start 
deteetor or other eleetrie fields downsfream of the start deteetor, the Mu atom will move from the 
C-foil to the MCP2 with uniform veloeity (vmu)- Its TOF from the C-foil to the MCP2 
determined from the TOF speetra is given by 

AtJ3r = fMu-to, (2.2) 


where tMu is the position (most probable value) of the Mu peak. Therefore, the Mu kinetie energy 
right after erossing the C-foil is given by 
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(2.3) 


where /mmu is the Mu mass and L the distanee between the C-foil and MCP2. It can be assumed that 
both the and the Mu leaving the C-foil have the same kinetic energy E'er = With this 

assumption, the energy loss in the C-foil (independent on the muon charge state when leaving 
the foil) is simply given by 


F'loss — ^kmod F'cp , 


(2.4) 


where Vmod is the HV at the moderator and Vc the HV at the C-foil. 

The energy loss in the ~2 /rg/cm^ C-foil for various energies are summarized in table 1 
and are plotted in figure 4. From the energy loss, a stopping power of 5 = 0.52(6) keVem^/jUg for a 
muon energy of 12 keV is obtained which has to be compared with the value of 0.70(1) keVcm^//rg 
in [12]. The uncertainty of the extracted stopping power and its deviation from the value in [12] 
is related to the uncertainty of the C-foil area density (0.2 /rg/cm^) originating from the non¬ 
uniformity of the foil and the uncertainty in the thickness from production to production. 

To check for the correctness of the assumed distance L = 563 mm used in Eq. (2.3) we calcu¬ 
late the TOF of the from the C-foil to the MCP2 (At‘'+‘^) using simple kinematic calculations and 
compare it with the TOF of /r+ peak determined directly from the measured TOF spectra (At)]'^®), 
e.g. figure 2. In order to calculate the /r+ TOF from the C-foil to the MCP2, we need to consider 
the various regions defined by the distances di as shown in figure 5. Since fhe TOF measurements 
have been accomplished without any HV at the RA conical lens, L3 einzel lens and no electric 
fields in the sample region the traveling from the C-foil till the MCP2 experience only regions 


- 5 - 




Figure 4. Energy loss in the 10 nm C-foil determined from the Mu peak in the TOE spectra for various 
incoming muon energy at the foil surface. 


of constant electric potential and two small regions in the start detector of constant electric field. 
The non-relativistic equations of motions for the various regions di are given by: 
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where m is the /r+ mass, E, the /i+ kinetic energy at the entrance of the /-region (^i = Ecf), q the 
charge of the particle, and Vc = —3.38 kV, Vp = —3.19 kV and Vg = —3.56 kV the various HVs 
applied at the grids of the start detector. The TOP /, in these various /-regions can be calculated 
using these simple relations: 
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Figure 5. (a) Schematic view of the TOF measurement, (b) Schematic of the start detector with various 
grids to deflect the electrons and defined regions of constant potential or constant electric held. Distances di 
(i = 1 — 6) are defined for on the beam axis. 


where v,- is the velocity when entering the region di, and a; the acceleration in the region di 
of constant electric field. The total TOF is eventually given by the sum 

= (2.13) 

/—I 

As can be seen by comparing the 5-th with the 6-th rows of table 1, there is a very good 
agreement between At‘^+‘^ and confirming fhe consisfency of our TOF analysis, including 

fhe correcfness of all disfances and fhe assumption fhaf Mu and exiting fhe C-foil have same 
kinefic energies Ecf- 
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Table 1. Moderator potential ymod. incoming muon energy at the C-foil surface E^+, measured and cal¬ 
culated TOP (for the peak maximum) of and Mu from the start detector to MCP2. The energy of the 
particle after the carbon foil E’er is determined from . Knowing the PcF, the energy loss Pioss in the 

carbon foil can be calculated using Eq. (2.4). The assumes the distances di to d(, given by mechanical 
construction and also assumes mono-energetic hitting the C-foil with energy given by the moderator and 
C-foil electric potentials, E^+ = eV^nod + 3.38 keV. 


Rmod (kV) 

8.5 

10.0 

12.0 

13.5 

15.0 

16.5 

18.0 

(keV) 

11.88 

13.38 

15.38 

16.88 

18.38 

19.88 

21.38 

Ecf (keV) 

10.84(6) 

12.37(7) 

14.49(9) 

16.11(10) 

17.67(12) 

19.17(14) 

20.78(15) 

Floss (keV) 

1.04(6) 

1.01(7) 

0.89(9) 

0.77(10) 

0.71(12) 

0.71(14) 

0.60(15) 

(ns) 

131.4(5) 

123.0(5) 

113.7(5) 

107.8(5) 

102.9(5) 

98.8(5) 

94.9(5) 

* (ns) 

152.7(5) 

139.8(5) 

126.4(5) 

118.6(5) 

112.1(5) 

106.8(5) 

101.6(5) 

(ns) 

152.9(3) 

140.0(3) 

126.5(3) 

118.4(3) 

112.0(3) 

106.7(3) 

101.8(3) 


3. Geant4 simulation of the energy losses in a thin carbon foil 

A Geant4-based (version 9.4 patch 04) simulation of the beam propagation in the LEM beam line 
has been accomplished. Figure 6 shows the geometry implemented in the Geant4 simulation. Beam 
line components relevant to this study are included, from the LI einzel lens, spin rotator (SR) to the 
sample chamber. The detailed geometry including radiation shield, sample holder and cold hnger 
of the cryostat are also implemented. 


Moderator 



Figure 6. Schematic view of the LEM beam line as implemented in the Geant4 simulation. The beam which 
is given by the black solid line is simulated starting after the moderator. The arrows indicate spin (blue) and 
momentum (red) directions. 

Precise electric and magnetic held maps are inputs to the Geant4 Monte Carlo simulation. The 
magnetic held map of the SR was measured in a volume 5x5x5 cm^ around the origin of the 






































































































SR coodinate system. The eleetrie field of the SR was ealeulated using the eommereial OPERA 
finite element programs (TOSCA/OPERA-3D) [26], and the eleetrie field maps were ealeulated 
with the finite element software COMSOE [27]. Due to the modified sample plate setup in this 
work deseribed in See. 4, the eleetrie field maps of the eonieal lens (RA) and the eopper sample 
plate with eylindrieal nose were re-ealeulated. The eleetrostatie module of COMSOE was used and 
a fine mesh was applied for higher aeeuraey ealeulations. A 2 mm spaeing of grid points are used 
for the eleetrie field maps. An example of the eleetrie potential map in the sample region is shown 
in figure 10. Initial eonditions of the jj.^ beam is summarized in table 2. The typieal number of 
events generated is ranging from 10^ to 10^, sueh that the statistieal uneertainty is redueed to less 
than 1%. 

The low-energy physies proeesses eurrently not available in Geant4 have been implemented 
to deseribe the energy losses and Mu formation proeesses in the thin C-foil of the start deteetor. 
When a /r+ is impinging on the C-foil our Geant4 simulation performs following operations: 

• Deeide about the eharge state of the exiting muon, between /r+. Mu and Mu^ using the 
yields given in figure 3. 

• Caleulate the energy loss (same for all partiele eharge state) assuming a Eandau distribu¬ 
tion with most probable value (MPV) given by the energy loss determined from the TOE 
measurements. 

• Compute angular seattering using the standard Geant4 paekage for multiple-scattering. 

More details regarding these three operations are given in the following subsections. 

Formation of muonium in the carbon foil 

Eormation of “foil” Mu at the thin C-foil is implemented by velocity scaling of existing data from 
proton-C-foil data [11, 13, 14] as shown in figure 3. These “foil” Mu will be stopped when they 
reach a material interface. 

Muon energy loss in the carbon foil 

In Geant4.9.4, models simulating the energy loss and its fluctuation are implemented in the 
C-i-i- class G4MuIonisation. By default, for energy below 200 keV, G4BraggModel is used 
where the energy losses are derived from the tabulated stopping power for proton using velocity 
scaling [15]. Energy loss fluctuation of is simulated by means of the G4IonFluctuation model. 
Eor a thin absorber, the energy fluctuation is based on a very simple two energy-level atom model 
and the particle-atom interaction give rise either to an atomic excitation or an atomic ionization 
with energy loss distributed according to [16]. 

In our simulation, the energy loss is simulated based on the values determined from the 
TOE measurements and its fluctuation is implemented using a Eandau random number generator 
based on CERNEIB [17]. It was found by E. Eandau that a certain linear function of the energy loss 
has, under certain assumptions, a universal (i.e. parameter free) density [18]. The generated random 
number from the universal Eandau distribution X is first shifted to have only positive energy losses 
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Figure 7. Energy loss distributions of 15.38 keV muons in a 10 nm thick C-Foil with a density ~2 /rg/cm^, 
obtained using the default Geant4 package (dotted black line) or our extension based on Landau distributed 
energy losses (red solid line) where the MPV energy loss (Eioss = 0.89) is taken from the TOE measurements. 
The energy loss spectrum extracted from the TOE spectrum is also shown (blue triangles). 


{X + 3.5) and then scaled linearly so that its MPV coincides with the measured ^loss in table 1, i.e. 
the randomly generated energy loss distribution is 

(3.1) 

Here, it is thus assumed that the energy loss distribution goes down to zero. 

A comparison between energy losses extracted from the TOP spectrum and simulated energy 
losses using standard or our extended Geant4 version are shown in figure 7. A better agreement 
between simulations and measurements is achieved when parameterizing the energy losses using 
the Landau distribution of Eq. (3.1). Interestingly, the Geant4 default simulation gives a better 
agreement below 0.5 keV. A cutoff approach was tried, however no improvement was achieved for 
the fitting of the TOE spectra. It is important to stress that the main attention is on the high-losses 
tail because it impacts the first few 100 ns of muSR measurements. 


Muon multiple scattering in the carbon foil 

In previous versions of Geant4, the measured transmissions of /r+ beam from the start detector till 
the sample were poorly reproduced due to the underestimation of the multiple-scattering process 
from C-foil [14]. However, recent versions of Geant4 have better physics models of multiple¬ 
scattering which reproduce correctly the Meyer scattering [13,20]. In this paper, multiple Coulomb 
scattering is simulated by using G4MuMultipleScattering based on the model from Urban [21]. 
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Figure 8. Measured and simulated TOR The green curve represents the contribution arising from “unmoder¬ 
ated” muons which has been added to the simulated TOR as an empirical sum of two “Lorentzian” to achieve 
good matching with the data (see text for details). 

Validation of the simulations 

To validate the implementation of the low-energy processes, the simulated TOR spectra are com¬ 
pared to the measured TOP spectra. In figure 8 a comparison between simulated and measured 
TOP is presented. The green curve corresponds to the contribution of these “unmoderated” muons 
which have to be assumed in order to match the measurements with the simulations. These are 
muons leaving the moderator not as epithermal muons at eV energy but as the non-fully moderated 
tail of the muon beam with keV energies. 

About 40% of the /rE4 beam hits the moderator target where about one half is stopped [5, 6]. 
This means that about 20% of the incoming beam will go through the moderator as “unmoder¬ 
ated” /!+ with a mean energy of several hundred keV and a low-energy tail ranging down to few 
keV energies [22]. Even though a large fraction of the “unmoderated” will not be reflected by 
the electrostatic mirror, they still contribute about 10-15%, depending on the moderator HV, to the 
total /f + which are impinging on the sample. It has been shown [23] that the electrostatic mirror 
which is set at the same HV as the moderator (Vnod) deflects particles with kinetic energies {E^) in 
the range eV^od < Ek < 2eVniod by 90°. 

Hence, “unmoderated” muons with slightly larger kinetic energy compared with the “mod¬ 
erated” muons are deflected towards the sample region. This explains why the “unmoderated” 
fraction contribute to the time spectrum at slightly earlier times compared with the “moderated” 
muons as well visible in Eig 8. 

The existence of “unmoderated” muons has been experimentally verified by taking data with¬ 
out any Ar-N 2 layer at the moderator. The corresponding TOE spectrum is shown in figure 9. 
Differently from the situation in figure 8 this measurement was performed with non-zero voltage at 
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Figure 9. (a) TOF spectra of the “unmoderated” determined from a measurement without any Ar-N 2 
layer on the moderator for various moderator and electrostatic mirror HVs (solid black line —12 kV, dotted 
red line —13.5 kV and dashed blue line —15 kV). (b) Measured “unmoderated” TOF at 12 kV beam transport 
(blue dotted line) fitted with the sum of two Lorentzian functions (red solid line). For comparison a TOF 
spectrum (black dashed line) of “moderated” muons measured with the Ar-N 2 moderator rescaled by a factor 
of 10 is also shown. 

the conical lens. The contribution of “unmoderated” muons to the TOF spectra can be empirically 
described by the sum of two Lorentzian functions which account for the Mu and the peaks. The 
“unmoderated” muons TOF were poorly reproduced by Geant4 simulation even after implementing 
the Landau energy loss distribution described in Sec. 3. This is due to the insufficient knowledge 
of the phase space of these muons after exiting the moderator. An extensive study is required to 
reproduce the shape of the TOF. 

The “unmoderated” muon contribution given by the green curve in figure 8 results from the 
sum of two “Lorentzian”. The relative widths and amplitudes of this two “Lorentzian” peaks can 
not be assumed from the “unmoderated” measurement because of the different HV settings of the 
conical lens focusing the beam on the MCP2 and therefore are free parameters. In conclusion, sim¬ 
ulated and measured TOF spectra agree very well together if a Landau distributed energy straggling 
in the C-foil is used and a small fraction of “unmoderated” muons is accounted for. It is important 
to note that the TOF spectra of figure 8 at times around t E [108; 112] ns and t E [118; 122] ns 
cannot be reproduced simply by modifying the muon energy losses in the C-foil assumed in the 
Geant4 simulation. The data in these two regions can be reproduced only by the contamination of 
“unmoderated” muons. 

Simulations can then be used to determine the and Mu kinetic energy distributions and 
related arrival time distributions at the sample plate for any beam line settings (moderator HV, 
sample HV, conical lens HV etc). Slow and Mu tails cause detrimental distortions of the 
measured /iSR time spectra which need to be accounted for when considering the “early” part of 
the measured time spectra. The starting point of the time window where /rSR fit can be reliably 
applied without being distorted is dictated by the low-energy tail. 

The knowledge of the fraction of “foil” Mu entering the LEM spectrometer is also an essential 
input for the analysis of the /iSR data. As the precession frequency of the muon spin in the Mu 
atoms is a factor of 100 larger than that for a free muon [24], the “foil” Mu produced at the C-foil 


- 12 - 














(b) Upstream detector Downstream detector 



Figure 10. (a) Electric potential contour plot in the “nose sample plate” region calculated with COMSOL 
multiphysics [27]. The “nose sample plate” was designed for the dedicated experiment of [25]. The normal 
sample plate for jiSR experiments is a flat disk, (b) The sample plate region is surrounded with upstream 
and downstream positron counters. 


give rise to depolarization effects and reduction of the observable total decay asymmetry. This is 
because the Mu eventually stop not only in the sample region with a well defined constant magnetic 
field buf also on fhe fhermal shield of fhe sample cryosfaf and ofher elemenfs which may have 
depolarizing effecfs and are placed af various B-field values. In fhe simulation, we assume fhaf 
once Mu has formed, depolarizafion occurs independenfly of ifs sfafes (single! or friplef). 

The recenf upgrade of fhe LEM beam line was characferized mainly by fhe insertion of fhe spin 
rofafor and fhe moving of fhe frigger defecfor closer fo fhe sample region. The closer placemen! 
of fhe trigger detector to the sample region has caused a larger fraction of the “foil” Mu to reach 
the sample region, enabling a better study of the process related with “foil” Mu production. The 
insertion of the spin rotator opened the way for longitudinal /iSR measurements, broadening the 
spectrum of possibilities available at the PSI-LEM spectrometer. In addition, the spin rotator was 
designed also to reduce beam contamination into the sample region. However the insertion of the 
spin rotator changed in an still not fully understandable way the beam propagation, degrading the 
beam size at the sample position. It is probably the insufficient knowledge of the fringe fields of fhe 
newly inserted spin rofafor which do no! allow an exacf simulation of fhe fransporf of fhe muons. 


4. Upstream-downstream asymmetries and the muon beam sizes 


The beam size af fhe sample posifion is needed fo normalize and analyze fhe /iSR dafa buf also fo 
undersfand fhe EEM beam line and validate fhe Geanl4 fransporf simulation. Eor example knowl¬ 
edge of fhe beam size af fhe sample posifion is used fo remove fhe confribufions arising from muons 
no! impinging on fhe sample of inferesf. The so-called upsfream-downsfream asymmefry Aud can 
be used fo infer fhe beam spof size af fhe sample position: 


^ud(t) 


Nu{t)-Nd{t) 
Nu{t)+Nd{t) ’ 


(4.1) 


where A^u.d(0 are fhe fofal number of decay position defected as a funcfion of time, in fhe upsfream 
and downstieam defectors surrounding fhe sample region as shown in figure 10(b). The values of 
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Figure 11. (a+c) The black dots are the simulated upstream-downstream decay asymmetry And for various 
RMS values Xrms of the muon beam spot. These plots are used to determine the experimental muon beam 
sizes at the sample position. In the simulation. Mu is generated at rest on top of (a) a 1 mm thick 20 x 20 mm^ 
fused quartz disk on top of the nose sample plate, for year 2011 setup, without spin rotator and (c) the 
Ag coated nose sample plate, year 2012 setup, with spin rotator in the low-energy muon beam line. The 
horizontal dotted lines are the measured decay asymmetries, (b-td) Muon beam spot at 14 keV implantation 
energy measured with the MCP2. The top (bottom) panels are before (after) the LEM upgrade. Quartz and 
Ag samples were used for the determination of muon beam sizes because Mu emission into vacuum is absent 
in these materials and Aud are time independent. 


Aud given in this paper were obtained by fitting Aud(0 with a constant function for times larger 
than 200 ns. The fitted Aud for various implantation energies are shown in figure 11(a) and (c). 
The upstream-downstream asymmetry and its time evolution is also the central ingredient of the 
longitudinal /iSR technique. 

A MCP plate at the sample position can be used to perform not only measurements of the 
muon TOP, but also measurements of the beam profiles which are shown in figure 11(b) and (d). 
The MCP measurements provide a 2-dimensional profile of the muon beam, but it can only be 
used when there is no HV applied to the sample plate. On the contrary, the upstream-downstream 
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asymmetry is strongly correlated with the beam size, and it can be measured on-line and for any 
sample-plate HV. This plays a crucial role especially at low energy, when a high positive HV has to 
be applied to the sample plate (Sample) to tune the /r+ implantation energy ^implant which is given 
by 

^implant — ^l^mod ^f^sample ^loss • (4-2) 

The electric potential ensued by the HV at the conductive sample plate of figure 10(a) shows a 
curvature of the equipotential lines which gives rise to a radial force causing a defocussing of the 
muon beam and thus an increase of the muon profile at the sample plate. Note that this effect is 
particularly relevant for the “nose sample plate” shown in figure 10 which was developed for a 
dedicated experimenf, which looked for fhermal Mu emission info vacuum from mesoporous silica 
fargefs [25]. The sfandard /rSR sample holder, being a simple plate, do nol show such a sfrong 
curvafure and Iherefore fhe beam defocussing effecf is smaller. 

A large variation of fhe beam size af fhe sample position for small variation of fhe beam 
parameters when using Ibis “nose sample plate” has been observed. Because of fhis sensifivify 
a sfudy of fhe Aud asymmefry using Ibis “nose” sample plafe for various beam line sellings was 
performed lo invesligale fhe validity of fhe Geanl4 beam Iransporl simulalion, which is also relevanl 
for fhe analysis of fhe experimenf for fhermal Mu emission info vacuum [25]. 

The black lines in figure 11(a) and (c) shows fhe correlation belween fhe beam size and fhe 
asymmefry Aud- They have been compuled assuming muon decaying from fhe sample plate wilh 
a given Iransverse spalial dislribulion described by a 2D Gaussian function wilh a widlh xrms = 
JRMS = Nole lhal Aud can be larger lhan fhe Iheorelical-maximum-decay asymmefry of 0.33 
due lo fhe shielding effecl of fhe nose sample plafe on fhe downslream defector, i.e. position has 
a lower probability of reaching fhe downslream deleclor (according to Eq. (4.1), Aud —> 1 when 
fVu>fVd). 

By comparing fhe measured asymmetiies (horizonlal dolled lines) wilh fhe asymmefry versus 
beam size predicled from fhe simulations, fhe beam size can be exfracled. On his torn fhis beam 
size can be compared wilh fhe beam size oblained from a Iransporl simulation of fhe full LEM beam 
line slarling from fhe moderator fill fhe sample region, including fhe processes in fhe C-foil and fhe 
eleclric fields in fhe modified sample region shown in figure 10(a). The asymmefry measuremenf 
can be Ihus used to validale fhe Geanl4 beam Iransporl of fhe LEM beam line. 

Eigure 11(a) and (c) show a decrease of fhe asymmefry for decreasing implanlalion energy 
revealing lhal fhe beam size increases considerably wilh decreasing energy. This has lo be related 
to a subslanlial defocussing effecl when Ihe HV al Ihe sample is increased due to slronger eleclric 
fields and slower muon velocity. 

Erom Ihe correlation line deduced from Geanl4 simulation as depicted in figure 11(a), Ihe 
muon beam size al 14 keV implanlalion energy has a RMS value = 6.9 mm. This value 
compares well wilh Ihe value measured wilh Ihe MCP2 of 6.3 mm (The MCP2 is placed 18.5 mm 
downslream of Ihe nose sample plate and hence a slighlly smaller beam spol is expected) and wilh 
Ihe value of 7.2 mm computed wilh Geanl4 beam Iransporl. 

However for Ihe 2012 measuremenls, after Ihe LEM upgrade, Ihe beam size values oblained 
from Ihe MCP measuremenls and Ihe values from Ihe asymmefry measuremenf do nol agree wilh 
Ihe beam size oblained from Ihe Geanl4 simulation of Ihe beam propagation in Ihe LEM beam line. 
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Table 2. LEM beam line initial settings for the Geant4 simulations and initial beam parameters. 



High voltage setting (kV) 

Moderator 

12 

15 

Einzel lens, El 

7.19 

8.99 

Mirror 

12 

15 

Spin rotator, 

2.09 

2.29 

Spin rotator, 

62.3 G 

69.7 G 

Einzel lens, L2 

8.38 

10.484 

Lens, L3 

8.99 

11.483 

Conical lens, RA 

9.15 

11.9 


Initial 1 

3 eam parameters 

Beam size, Gx,y 

7.5 mm 

Beam divergence, ay y 

2 .0° 


This is calling for a verification and detailed investigation of the LEM beam line simulations. To 
study the sensitivity of Aud to various beam line settings the sample holder of figure 10 is used. The 
beam line settings of table 2 are for the SR setup, experimentally optimized. 

The dependence of Aud on the various beam line parameters has been investigated and sum¬ 
marized in figure 12. 

• (a) Muon implantation energy: Aud is increasing with decreasing implantation energy due 
to defocussing effects caused by the electric potential applied at the sample holder (see fig¬ 
ure 10(a)). The musrSim (new) represents the Geant4 simulation where the muon energy 
losses at the C-foil are parametrized using Landau distributions. Lrom figure 13, it can be 
seen that for a lower implantation energy, there is a considerable amount of which do not 
hit the “nose sample plate” which has a radius of 15 mm. Thus at lower implantation energy, 
there is a higher probability of decay positron being detected by the downstream detector, 
resulting in a reduced Aud asymmetry. 

• (b-i-c) Initial phase space of the muon heam: Initial phase space {xx',yy') after the accel¬ 
eration section of the moderator is taken to be Ox = Oy = 7.5 mm for the beam size and 
Ox' = Oy! = 2.0° for the beam divergence from a recent simulation [13]. The initial polar¬ 
ization vector is chosen as = (0.9848,0,0.17365) since the spin is rotated by 10° 
clockwise after traversing the electrostatic separator of the /xE4 beam line before the /f + is 
focused on moderator (see figure 1). 

• (d) Electric and magnetic field of the spin rotator: The magnetic field of the spin rotator 
was fixed to = —62.3(—69.7) G for 12(15) keV transport energy to obtain the experi¬ 
mental /x+ spin rotation of 20° counter-clockwise. The spin angle with respect to the jj.E4 
beam direction (corresponding to -x direction in the simulation) is then changed from -i-10° 
(the angle after the separator in the /rE4 beam line) to -10°. The electric field was varied 
from = 1.85 to 2.25 kV/mm and is chosen such that the beam is centered on the original 
sample plate position (16 mm downstream of the nose sample plate) as shown in figure 14(a). 
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Figure 12. Measured and simulated Aud versus /r+ implantation energy. Gxy and a^/yi are the standard 
deviation of the muon phase space after the acceleration section of the moderator, the electric field in 
the spin rotator and Vra the HV at the conical lens. For each plot only one parameter is varied whereas the 
other one are given in table 2. 


From figure 12(d) it is visible that that Aud is maximal at around 2.09 kV/mm. This occurs 
when the beam is centered on the sample plate (see figure 14(a)) and the downstream detec¬ 
tors are shielded from the positron by the material of the “nose sample plate” itself. 

The design value of 2.68 kV/mm which was expected from simply considering the relation 
V = ^ has thus to be tuned to 2.09 kV/mm to center the beam on the axis. This indicates that 
there are not well understood imperfections of the electric and magnetic fields that cause the 
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Figure 13. Simulated energy of the muon beam at sample-plate plane as a function of distance R from 
the beam axis for 2 different muon implantation energies: (a) 11 keV and (b) 3 keV. Region 1 is mainly 
contributed by muons with kinetic energy above 9.15 keV before reaching the RA. Region 2 is mainly given 
by muons with energy around 9.15 keV (= Vra) before reaching the RA. Region 3 is originated from Mu, 
which are not decelerated by the electric helds in TD and are not affected by the RA and sample plate electric 
potentials, and hence have higher energies than the /i+. Since Mu are not focused by the RA lens they have 
a distribution with larger transverse extension. 


/!+ to deviate from its original path. 

• (e) Electric potential of the conical lens (RA): If the focusing power is not optimal the Aud 
decreases since the beam becomes larger. 

• (f) Z-offset of the sample plate position: It is obvious that there is a strong dependence of 
the Aud on the z-position of the “nose sample plate”. Moving the sample plate downstream, 
in z-direction, will decrease the Aud as the downstream detector will be less shielded by the 
sample plate from positrons. 

Summarizing, the new musrSim with Landau distributed losses give asymmetries closer to 
the measured one. However, the agreement is not yet satisfactory. Some parameters could be 
slightly tuned around the design value to decrease the deviation between measured and simulated 
Aud asymmetries. The newly inserted spin rotator optimal settings do not correspond to the design 
value manifesting that some uncontrolled beam distortion and beam transmission may occur at 
this beam line element. Contribution of “unmoderated” /i+ is shown in figure 14(b). As they 
have a larger beam spot {xrms and yRMs) as shown in table 3 and hence a lower value of Aud- 
The larger beam spot could be caused by chromatic aberration in the beam transport optics due 
to the wide energy distribution of the “unmoderated” /r+. But this is not enough to explain the 
observed discrepancy between experimental and simulated Aud asymmetries as they accounted for 
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Figure 14. (a) Simulated beam mean position {x,y) at the sample plate as a function of the electric field in 
the spin rotator The F™ value in table 2 is chosen such that the beam is centered on the original sample 
plate position (16 mm downstream of the nose sample plate), (b) Measured Aud versus implantation 
energy with and without Ar-N 2 layer at the moderator. “Unmoderated” muons at the plate position have a 
larger beam size. 


only 10-15% of the measured time spectra. However, as shown in figure 12, by tweaking different 
parameters one could get agreement. We refrain at this stage to perform a multi-parameter tuning 
because of the complexity and the correlation between the various parameter. 

Table 3. Measured /r+ beam spot xrms and jrms for the “unmoderated” /r+ and “moderated” /r+ using the 
MCP2. It should be noted that due to the hnite size of the active region of the MCP2 (42 mm in diameter) 
the actual xrms and jrms could be larger, especially for the “unmoderated” /r+. 


Energy (keV) 

Spin rofafor 

Moderator 

Xrms (mm) 

jRMs (mm) 

fype 

15 

No 

No 

6.1 

6.0 

unmoderafed 

15 

No 

Yes 

5.3 

4.9 

moderated 

15 

Yes 

No 

7.2 

7.2 

unmoderafed 

15 

Yes 

Yes 

6.2 

6.0 

moderated 


The reason why there was a good agreement between the simulated and measured Aud prior 
to the LEM beam line upgrade, is related with the smaller beam size at the sample plate which 
imply a reduced defocussing effect at the “nose sample plate”. The insertion of the spin rotator has 
caused a degradation of the beam quality and an unexpected increase of the beam size at the sample 
position. This increase could be even more substantial for the “unmoderated” muon component. 
The larger RMS values (« 1.2 mm) for “unmoderated” muons could be sufficient to explain the 
observed smaller asymmetries in figure 12. 

5. Conclusions 

Simulations of fhe complefe LEM beam line afler fhe 2012 upgrade have been presenfed. TOE 
measuremenfs have been used fo calibrafe fhe energy losses in fhe sfarf defector. Excellenf agree- 
menf befween fhe measured TOE specfra and simulations has been reached only by using Landau 
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distributed energy straggling, aeeounting for muonium produetion, and aeeounting for a eontribu- 
tion of “unmoderated” muons with slightly larger kinetie energy whieh are parasitieally transported 
by the beam line. The measured stopping power in the C-foil eompares well with previous deter¬ 
mination. This good agreement between measured and simulated TOP speetra implies also the 
eorreetness of the assumed muonium yield in the C-foil whieh has been implemented using veloe- 
ity sealing of proton data. 

Detailed analysis of the /xSR time speetra require information of the muon arrival times, muon 
implantation energy and depolarization effeets related with energy losses and Mu produetion in the 
C-foil. These information ean be determined now for all beam line settings and muon implantation 
energies using the new Geant4 simulations. 

The beam spot size at the sample position is also a very important parameter when analyzing 
the /rSR data. Normalizations and total measurable deeay asymmetries depend on this parameter. 
This is even more important for longitudinal /rSR whieh is now possible due to the insertion of the 
spin rotator. Therefore the beam spot size at the sample position has been investigated by means of 
the Aud asymmetry using a dedieated sample plate (“nose sample plate”) with inereased sensitivity 
to beam ehanges. These studies have revealed a problem with the beam transport in the LEM beam 
line related with the newly inserted spin rotator and the parasitie transport of “unmoderated” muons 
whieh has eahed for detailed studies and hardware improvements. 
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